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Abstract 



The present paper reviews recent achievements 
O^on the ab initio determination of effective model 
^Hamiltonians aimed at the description of strongly 
CNcorrelated materials. These models (Heisenberg, 
xpf-— J, extended Hubbard, Kondo, etc) are crucial 
OtP solid state physicists for the description and un- 
"+^der standing of the electronic remarkable proper- 
ties such as magnetic orders, photo-induced mag- 
netism, transport properties, high t c super con- 
^ductivity, etc. We will see how the association 
Qatid the control of embedding techniques, fragment 
^spectroscopy, effective or intermediate Hamilto- 
>nian theories, provides a systematic and well con- 
strolled method for the determination of trustwor- 
thy models Hamiltonians. 

a ' 

_ _ i 

1. Introduction 



Since the synthesis of the first organic con- 
ductors in the late 70s with the TCNQ-based 
compounds^ and later the Bechgaard saltsS, the 
chemists have synthesized a large number of ma- 
terials presenting new and remarkable electronic 
properties. Besides the organic and organo- 
metallic metals, to cite only a few among the 
most attractive systems of the last decades, one 
can remember the spin chains^ and ladders^, the 
cnnner-rmHe hifh T„ siiner-conHiictors^. manfan- 



ites with giant magneto-resistance^, photo-induced 
magnets' 7 , etc... These materials span a very large 
variety of structural arrangements, phases and 
physical properties. Indeed one can find among 
them molecular crystals, transition metals oxides, 
Prussian blue analogs, etc. They can be isotropic 
as well as present large anisotropies both in their 
properties and structural arrangement, from the 
cubic, 3D, isotropic photo-induced magnets, to the 
2D, layered copper oxides or manganites, up to 
the quasi-lD stacks of the Bechgaard salts or ox- 
ides spin chains and ladders. These compounds 
present a wide variety of phases going from metals 
and super-conductors to semi-conductors, Peierls 
or spin-Peierls insulators, Mott insulators, charge 
or spin density wave insulators, they may present 
long range order such as anti-ferromagnetic order- 
ing, have exponentially vanishing correlation func- 
tions, etc. 

In view of this very large spectra of chemical 
compositions, structural arrangements and physi- 
cal properties one may wonder what all these com- 
pounds have in common. In fact all these ma- 
terials have in common an essential characteris- 
tic : the electrons responsible for their low en- 
ergy physics are few (per unit cell) and local- 
ized both spatially and energetically. The spa- 
tial localisation of these, from now on called, 
active electrons induces a strong correlation be- 
tween their relative nositions ( r>n tr-i . fo 
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pi (fi ; rJJ pi (f 2 ; f^)-pi (fi ; r* 2 ) Pi (^2 ; r\), where 
P2 is the two electrons density matrix and pi is the 
one electron density matrix) and the impossibility 
to describe their movement — even in a qualita- 
tive manner or at the zeroth order — by an inde- 
pendent electrons representation. These systems 
belong to the strongly correlated universality class 
and essentially mono-electronic methods such as 
tight-binding, mean-field approximations or even 
density functional methods fail to describe their 
electronic structure^. The active electrons are en- 
ergetically localized close to the Fermi level and 
are thus responsible for the low energy properties 
of the materials. It is thus convenient to parti- 
tion the electrons of these strongly correlated ma- 
terials in a similar way as the usual core-valence 
separation in atoms. On can define 2 classes, the 
core electrons, essentially localized on the energeti- 
cally deep, closed-shell orbitals of the basic entities 
(atoms or molecules according to the chemical na- 
ture of the crystal), and the previously cited active 
electrons supported by localized, partially-filled or- 
bitals, presenting charge and/or spin fluctuations 
in the ground and low energy-excited states. 

While in finite, small systems, the treatment at 
the ab-initio level of strong electron-correlation ef- 
fects can be achieved in a satisfactory manner, in 
infinite systems this is still an essentially unsolved 
question. The difficulty is even larger when ex- 
cited states are involved or thermodynamic prop- 
erties seeked. In order to by-pass this problem and 
to understand the microscopic mechanisms under- 
lying the spectacular properties of these materi- 
als, physicists have developed the usage of model 
Hamiltonians on the same bases than the semi- 
empirical models used by quantum chemists in the 
50's and 60's^ — Zero Differential Overlap (ZDO) 
approximation^ and reduction to the dominant in- 
teractions — but with somewhat simpler expres- 
sions. Among the most famous models one can cite 
the spin models such as the Heisenberg-Dirac-Van 
Vleck model o 11 ! 12 for magnetic half-filled systems, 
the t — J models^ for doped magnetic systems, the 
Hubbard^ and extended Hubbard models — sim- 
plified versions of the Pariser-Parr-Pople^ approxi- 
mation — that can describe metal-insulator phase 
transitions. A simple description of these models 
is given in the appendix. 

The description and understanding of the elec- 
tronic structure, low energy physics and the asso- 
ciated macroscopic properties of these materials is 
therefore achieved within a three steps procedure. 



• The first step is the determination of an ap- 
propriate model. Such a model, aimed at de- 
scribing the low energy physics of the mate- 
rial, must define the active electrons, the ac- 
tive orbitals and the form of the dominant 
interactions between them. 

• The second step consists in the evaluation 
of the amplitude of these interactions. This 
evaluation can, in some well conditioned 
cases, be extracted from experimental data, 
however in most cases there is no unique and 
unambiguous way to determine such param- 
eters from spectroscopy or thermodynamic 
properties measurements. 

• The last step is devoted to the determina- 
tion of the collective properties of an assem- 
bly of electrons on a lattice, described by 
the previously-defined, effective, local inter- 
actions . 

We will see in this paper how it is possible to 
construct, step by step, a coherent procedure for 
the determination of trustworthy models describ- 
ing the low energy physics of strongly correlated 
materials. We will also see how quantum chem- 
istry can play a major role, both in the deter- 
mination and validation of a model and in reli- 
able evaluations of the effective interactions ampli- 
tudes. The last step is usually devoted to physi- 
cists since it requires the usage of specific ana- 
lytical methods such as renormalisation theories- 
bosonisation 15 , Bethe ansatz 16 , etc as well as spe- 
cific numerical methods such as Density Matrix 
Renormalisation Group (DMRG) methods 17 , ex- 
act diagonalisations, Monte-Carlo, etc. It aims at 
the determination of the macroscopic order from 
the microscopic effective model. It is however fas- 
cinating to follow the whole process and under- 
stand the intrinsic link between the local chemical 
structure of the materials and their macroscopic 
properties. 

2. Requirements for a good 
model 

Before any discussion on how to built an accu- 
rate and reliable model for strongly correlated ma- 
terials one should examine what means good, accu- 
rate and reliable, that is what do we expect from 
such a model. The model is supposed to fulfill 
several properties that can be ordered into three 
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sets : the formal properties, the physics we would 
like to reproduce and the transferability of the in- 
teractions. 

2-A. Formal properties 

These requirements are quite obvious and have 
already been mentioned. They can be summarized 

as 

• Hermiticity of the model Hamiltonian. 

• Small number of electrons and basis set or- 
bitals per unit cell. 

• Small number of basic interactions. 

2-B. Modeling the low energy physics 

The final aim of the model is to be able to 
produce a reliable representation of the low en- 
ergy physics and the properties of the material, 
that is to be able to described and interpret not 
only phase transitions but also experimental re- 
sults such as low energy spectroscopy, low en- 
ergy photo-emission or photo-absorption, NMR, 
X-ray or neutron diffraction, magnetic susceptibil- 
ity, conductivity, etc. We therefore need not only 
an accurate reproduction of the first excitation en- 
ergies but also accurate reproduction of the asso- 
ciated wave-functions. In particular it is crucial to 
well reproduce the charge and spin arrangements 
as well as the relative weight of the major con- 
figurations in the ground and low energy excited 
states. This means for instance, to reproduce ac- 
curately the ratio between the coefficients of the 
Valence-Bond (VB) neutral and ionic configura- 
tions, of the different spin arrangements within the 
neutral configurations, etc. Such requirements on 
the wave function suppose that 

• the low energy states we would like to re- 
produce (let us call them the target states 
l^i), i = l---m t ) are clearly defined (for 
instance states essentially based on neutral 
Valence Bond configurations) and based on 
a small set of common active orbitals, 

• the projection of the exact target states wave 
functions onto the vector space supporting 
the model Hamiltonian (called the active 
space from now on, P being its projector op- 
erator) should be large (one should notice 
that there is no requirement for the active 
space to be complete) 



• the exact target states excitation energies 
(Ei, i = 1 ■ ■ ■ m t ) should be well reproduced 
by the model hamiltonian eigenvalues (e^) or 
a convenient subset of them, 

• the model Hamiltonian eigenstates (noted 

i = l---m a , m a > m t ), or a subset 
of them should mimic the projection of the 
target states onto the active space, that is in 
the ideal situation 

Vz = l---m t , |6> = P\^i)/^(%\P\%) 
2-C. Locality of the interactions 



As already mentioned, the model Hamiltonians 
are defined on a basis set of local active orbitals 
and thus the microscopic interactions are local, in- 
volve only a small number of electrons and a small 
number of centers, usually one or two. Let us cite 
for instance, the one-electron, two-centers hopping 
integral (as in the t — J or extended Hubbard 
models) or the two-electrons, two-centers super- 
exchange integral (as in the t — J or Heisenberg 
Hamiltonians). One should however point out that 
the concepts of center or site are not necessar- 
ily synonymous of atom, but rather the expres- 
sion of the basic building block participating to 
the low energy physics of the system. For supra- 
conducting copper oxides for instance, a center will 
be a copper atom, the oxygens and other atoms 
being excluded from an explicit representation in 
the model. For molecular crystal a center will be 
the whole building molecule (such as TMTSF— , 
TCNQ 37 , etc) and its HOMO and/or LUMO will 
be the localized active orbitals supporting the low 
energy physics of the materials. The concept of 
locality can thus refer to fairly large units from 
the molecular point of view but of size always in- 
ferior to the unit cell and therefore localized from 
the crystal point of view. The locality assump- 
tion has the immediate and important consequence 
that the interactions do not depend of the number 
of explicit active centers and that, provide a set of 
orbitals and associated energies, a n-centers inter- 
action will act in a finite system in the same way as 
in an infinite one. This property is usually called 
transferability. The hopping or the exchange in- 
tegrals can, for instance, be defined from the low 
excitation energies of a two centers system. The 
exchange integral is the singlet triplet energy dif- 
ference in a two active electrons system and the 
hopping integral is half the symmetric versus an- 
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tisymmetric doublet energy difference, in a single 
active particle system. 

The exchange integral J ab 

2 e - 2 orb. - 2 sites 

T I 
• • 

a b 
Jab = E(sg) - E(tp) 

Hopping integral t a b 

1 e - 2 orb. - 2 sites 

T — 
a b 

t ah = {E{db + )-E{db_))/2 

Despite their local definition in the model Hamil- 
tonians scheme, it is clear that these interactions 
are effective and one should wonder whether they 
are truly local or whether they do involve in an ef- 
fective manner non local mechanisms. This is the 
question behind the transferability problem. Its 
answer is crucial for the determination of reliable 
models. Indeed, without the assurance of a true 
locality of the interactions the whole model hamil- 
tonian scheme fails. The derivation of the macro- 
scopic properties of an infinite lattice is based on 
the idea that this infinite system is the limit of 
increasing large size finite systems. The implicit 
hypotheses behind it are that the interactions in- 
volve a finite and small number of bodies and that 
their nature and amplitude do not change when 
the system size increases, i.e. they are perfectly 
transferable. The effective inclusion of non local 
mechanism is incompatible with these hypotheses 
since it would render the effective integrals depen- 
dent on the system size. 

One should however remember that the exact 
Hamiltonian can always be expressed in an atomic 
basis set where all interactions are local and trans- 
ferable. The non locality of a model Hamiltonian 
can therefore always be fixed by an appropriate 
modification of the model, usually consisting in 
the addition of longer range interactions, interac- 
tions involving a larger number of centers or an 
enlarged definition of the active space. The trans- 
ferability property is therefore an useful tool for 
the verification of the model pertinence. A model 
which is not transferable usually neglects the ex- 
plicit treatment of an important physical effect and 
the transferability is restored by the inclusion of 



the omited mechanism within the model Hamilto- 
nian. The whole process is however limited by the 
complexity of the final model, the number of basic 
interactions and the number of bodies involved in 
it. 

3. Locality of effective 
interactions 

In this section we will investigate the locality 
problem for the major types of effective interac- 
tions found in model Hamiltonians, namely hop- 
ping, exchange and repulsion. In this analysis we 
will suppose that the orbital space has been parti- 
tioned into the 

core orbitals , these orbitals are essentially 
doubly-occupied in the target states and sup- 
port the inactive electrons, explicit reference 
to them is usually omitted in the model 
Hamiltonian, 

active orbitals , they support the active elec- 
trons responsible for the low energy physics, 
these orbitals are supporting the model 
Hamiltonian and have fluctuating occupation 
and/or spins in the target states, 

virtual orbitals , these orbitals are essentially 
empty in the target states and omitted in the 
model. 

The above orbital basis set will be supposed to be 
local, orthogonal and the orbitals composition (as 
can be given by their expansion on the atomic or- 
bitals basis set) and energies already defined. We 
will come back later on the problem of these or- 
bital and orbital energies definitions and analyze 
the non- locality contained in it. This partition of 
the all electrons basis set is crucial and will be ref- 
ered at all over this paper. 

The effective integrals can be decomposed (us- 
ing the above partition of the orbitals basis set) 
into a zeroth-order part corresponding to the ex- 
act Hamiltonian integrals within the active space, 
and a higher order part that takes effectively into 
account the effect of the configurations out of the 
active space. It is in this higher order part that non 
local mechanisms may be introduced through long- 
range dynamical polarisation or correlation effects. 

We will start by the analysis of the effective bi- 
electronic integrals since their case is somewhat 
simpler than for the mono-electronic ones. 
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3-A. The effective exchange integrals E 



The effective exchange integrals appear in mag- A 
netic models such as the Heisenberg or t— J Hamil- 
tonians. These models are suited for the descrip- 
tion of very strongly correlated systems where the 
on-site bi-electronic repulsion is so high that the A i 
explicit reference to configurations where an active 
orbital is doubly occupied can be excluded from 
the model. The essential role of the effective ex- 
change integrals is therefore to take into account 
the effects of these double occupancies on the inter- 
actions between neutral Valence-Bond configura- 
tions, that is to position properly the local singlet 
to triplet excitation energy between two nearby 
centers. 

Within the complete configurations space de- 
fined from the active orbitals (CAS), Anderson^ 
has clearly detailed this mechanism and provided 
a perturbative expression of the effective super- 
exchange integral J ab (see figured)) 
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FIG. 1: Anderson super-exchange mechanism 
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where t a b is the hopping between the magnetic or- 
bitals a and b and U is the Coulomb repulsion in- 
tegral when a magnetic orbital is doubly-occupied. 

Things are a little more complicated when the 
exchange mechanism is not direct but through a 
chemical bridge. This is for instance the case 
in most transition metal oxides where the super- 
exchange interaction between two d orbitals of 
nearby metal atoms is mediated through the oxy- 
gen atom(s). In this case the leading perturbative 
term comes at the fourth order and two different 
mechanisms add (see figure |2J). 
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FIG. 2: Through-bridge super-exchange mechanism. In this 
example we supposed that the process goes through the occupied 
orbitals of the bridge, represented by x. It may also go through 
the unoccupied orbitals or both. 
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where 



L bx 



l bx 

±ax 
l ax 

A x 
U 
A 2 



(xx\H\xb) is the hopping between orbitals 
b and x with a spectator electron on x 
(bx\H\bb) is the hopping between orbitals 
b and x with a spectator electron on b 
(axx\H\axa) is the hopping between orbitals 
a and x with two spectators electron on a and x 
E (|ax 2 6)) - E (\a 2 xbf 

E(\ax 2 b)) -E(\a 2 x 2 )^ 

E (\ax 2 b)) - E (\a 2 b 2 )] 



It should be noted that the second path 
goes through a configuration where the doubly- 
occupied orbital of the bridge has been totally 
emptied (or totally doubly-occupied if the pro- 
cess goes through the unoccupied orbitals of the 
bridge). Such a configuration is usually very en- 
ergetic and its contribution negligible in front of 
the first path. In that case the through-bridge su- 
per exchange integral can again be written in an 
Anderson's form J ab = kt 2 ab /U provided that an 
effective, through-bridge, hopping integral t a b be- 
tween the two magnetic orbitals is defined as (see 
figure El) 
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It is clear that all these processes included in the 
effective exchange integral are totally local. Dy- 
namical polarisation or correlation processes that 
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FIG. 3: Through-bridge hopping mechanism. In this example 
we supposed that the process goes through the occupied orbitals 
of the bridge, represented by x. It may also go through the 
unoccupied orbitals or both. 

may act on top of the them (for a review of 
all diagrams contributing to the unbridged super- 
exchange mechanisms at the second order of per- 
turbation see ref.— ) remain local as long as they 
involve only orbitals located on the magnetic or 
neighboring atoms, including the possible bridge. 
The only possibility to include non-local contribu- 
tions is therefore the dynamical polarisation and 
correlation processes where the excitations take 
place on a center far from the magnetic orbitals. 
The leading contributions of such terms (second 
order of perturbation) correspond to the super- 
exchange plus polarisation diagrams in the DeLoth 
denomination^ (see figure Ej) 
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FIG. 4: Leading non local contributions to the exchange effective 
integral, i and i* are respectively occupied and vacant orbitals 
on a distant site. 



The corresponding perturbative contribution is 

(ib\—\i*a)(i*a\Mib) 

' 1 ri2 1 ' ' 1 rig 1 ' 

AE 

which varies as a function of the distance R be- 
tween the magnetic distribution ab and the distant 
local excitation dipole-dipole interaction, 

that is as 1/ (R 3 ) 2 = 1/R 6 . 

One can therefore reasonably consider that these 
non-local contributions are negligible and that the 
exchange effective integral is really local. This lo- 
cality property has been numerically verified by 
Bias and coworkers 20 from ab-initio calculations 
on a variety of compounds such as the KMF3 
and K 2 MF 4 (M = Cu, Ni) family or hight T c per- 
ovskites. 



3-B. The effective coulomb repulsion integrals 



The effective coulomb repulsion terms act on the 
diagonal of the effective Hamiltonian. They corre- 
spond to screened bi-electronic repulsions. This is 
through this screening that non local effects may 
occur. 

The screening that may come from the active 
electrons is treated explicitly in the model Hamil- 
tonian, for instance through hybridisation and de- 
localisation processes, and therefore do not con- 
tribute to the effective repulsion integral. It re- 
mains the dynamical polarisation and correlation 
effects. 

The main non-local contributions come from 
dynamical correlation effects through double- 
excitations from two core orbitals towards two vir- 
tual ones. However, these excitations yield equal 
contributions (at the second order of perturbation) 
to all VB configurations in the active space and 
give only a shift in the definition of the model 
Hamiltonian zero of energy. Thus they are non 
pertinent and can be omitted. 

The remaining leading non-local contributions 
come from distant single-excitations. For the sake 
of simplicity, we will suppose that the core (i) 
and virtual (i*) orbitals are eigenvectors of the 
core Fock operator (as defined from the antisym- 
metrized product of all doubly-occupied core or- 
bitals). The leading term comes therefore from the 
coupling between the (i i*) distant distribution and 
the local (a a) distribution (if the two electrons are 
on the same orbital), (ab) or (ab) distributions (if 
the two electrons are on different orbitals). See 
figure EJ) • 



a a 




a a 



FIG. 5: Leading non local contribution to the on-site effective 
repulsion integral between two electrons in the same orbital. 
i and i* are respectively occupied and vacant orbitals of the 
distant excitation 



In the case where the effective repulsion integral 
acts between two electrons in the same orbital, the 
leading non-local contributions are at the second 
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order contribution 

(ia\— \i*a)(i*a\ — \ia) 

x 1 r 12 I ' \ 1 ri2 1 / 

This terms vary as a function of the distance R 
between the aa distribution and the distant exci- 
tations charge to dipole interaction, that 
is as 1/ {R 2 f = 1/i? 4 . 

When the repulsion integral acts between two 
electrons in different orbitals, this dominant non- 
local contributions vary as a dipole-dipole interac- 
tions, that is as 1/ {R 3 f = 1/i? 6 . 

The non-local contributions to the repulsion in- 
tegrals vary therefore either as 1/R 4 or as 1/R 6 , 
that is fast enough to consider that these effective 
integrals are essentially local. 

3-C. The effective hopping integrals 

The case of hopping integrals or other mono- 
electronic integrals (such as orbital energies for 
instance) is somewhat more complicated. Indeed 
in their zeroth-order description (projected exact 
Hamiltonian in the active space) these integrals 
are coupled to the Fock operator, that is both to 
the nuclear and the electronic charge repartition 
on the entire system. It is clear that, through this 
coupling, the mono-electronic integrals are non lo- 
cal. However, once provided the definition of the 
orbitals (core, active and virtual), their composi- 
tion in terms of the atomic orbitals and their en- 
ergies, the other contributions to the hopping in- 
tegrals can be analyzed as was done for the super- 
exchange or repulsion integrals. Let us decom- 
pose the problem into three terms, the first term 
is the coupling with the core orbitals (and the cor- 
responding charge density of the entire system) 
through the Fock operator, the second term is con- 
tributions from the active orbitals, the third term 
is the excitations out of the active space that are 
included in an effective manner. 

The first term is totally defined by the Fock op- 
erator. The problem intrinsically depends upon 
the average charge repartition on the entire system 
and is solved by the definition of the orbitals them- 
selves through a self-consistent, mean- field approx- 
imation. This is the major non local contribution 
and its consequences is that the infinite system is 
never really treated as the limit of true large but 
finite system — since the edge effects should not 
be taken into account in the orbitals and orbital 
energies definitions. A proper definition of the or- 
bitals therefore suppose that the finite systems are 



in fact embedded in a mean field approximation of 
the rest of the crystal. 

The second term reminds us that the hopping 
integrals are dependent on the specific configura- 
tions they act upon. For instance the transfer of 
an electron on a bond depend on the number and 
the nature of the spectators (active) electrons on 
it. 

1L _ _t_ _L 

a b tj b a b 

FIG. 6: hopping with a spectator electron in a. 



tab = (aa\H\ab) 

= (a\H\b) + (aa\l/ri2\ab) 

11_ L L A 

a b tjj a b 

FIG. 7: t°t hopping with two spectator electrons in a and b. 



= (aba\H\abb) 

= {a\H\b} + (aa\l/r 12 \ab) + (ba\l/r 12 \bb) 

The hopping integrals therefore depend upon the 
distant active space charge distribution of the con- 
figurations they act on. The corresponding contri- 
butions vary as (xa\l/r\2 — xb), that is as 1/R 2 
if R is the distance between the (ab) and the dis- 
tant (xx) distributions. This dependence should 
be explicitly include in the model hamiltonian def- 
inition. However, despite its rather slow decrease 
with R, it is usually ignored, except may be for 
the dependence to the occupation of the orbitals 
localized on the same sites as a and b. 

The third term contributions are included in 
an effective manner. Their leading non-local pro- 
cesses come from distant single-excitations acting 
on the transfer, as shown in figure |H1 
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i* ; 
i 

FIG. 8: Leading non local dynamical contribution to the effec- 
tive hopping t a b- 

The Fock operator being defined as previously, 
the perturbative expression of the dominant is 

(ai\l/r 12 \ai*)(ai*\l/r 12 \bi) 



(ai\l/r 12 \bi*)(bi*\l/r 12 \bi) 

which varies as 1/R 2 x 1/R 3 = 1/R 5 , that is fast 
enough to be considered as a non-significant source 
of effective non-local processes. 

As a summary the major types of interactions 
encountered in model Hamiltonians are essentially 
local in nature, provided that the coupling to 
the average charge repartition of the entire crys- 
tal (which is the major non-local mechanisms) is 
treated through a proper definition of the orbitals 
and orbital energies. It is therefore pertinent to 
evaluate the interactions amplitudes and check a 
model validity by ab-initio calculations on properly 
embedded fragments. Such a technique presents a 
lot of advantages. 

• All short range effects can be explicit and ac- 
curate treated by the inclusion of the active 
sites local environment in the computed frag- 
ment. 

• The dynamical processes that strongly affect 
the values of the effective interactions can be 
accurately evaluated. Indeed, it is possible 
to use large atomic basis set and up-to-date 
quantum chemical methods for the fragment 
spectroscopy. 

• One has access not only to the local excitation 
energies but also to the associated wave func- 
tions. In particular the knowledge of the rel- 
ative contributions of the major active space 
VB configurations is determinant for the good 
reproduction of the material properties by the 
model. 

The question of the embedding is therefore crucial, 



as it determines the composition of the orbitals 
(core, active and virtuals) as well as their energies. 



4. How to built a good bath 

According to the previous analysis a good embed- 
ding should be such that the fragment orbitals and 
their energies are similar to what they would be in 
the crystal. Reaching this goal necessitates that 
the effects of the rest of the crystal on the frag- 
ment are properly reproduced. These effects may 
be partitioned into long range effects and short 
range effects. 

The previous analysis showed us that the ma- 
jor long range effects are mono-electronic, comes 
through the Fock operator definition and more 
specifically through the coupling to the entire sys- 
tem average charge repartition. Considering that 
these strongly correlated systems are most of the 
time ionic, this term is large and corresponds to 
the Madelung potential. It can be easily repro- 
duced by a set of positive and negative charges re- 
spectively positioned at the location of the cations 
and anions of the rest of the crystal. The volume 
involved should be large enough and centered on 
the fragment. Border charges can be adapted ac- 
cording to the Evjen procedure in order to increase 
the potential convergence. 

The short range effects of the rest of the crystal 
on the fragment are essentially exclusion effects. 
They come from the orthogonality requirement of 
the fragment orbitals to the orbitals of the rest of 
the crystal. This effect is crucial since in the vac- 
uum the orbitals of the fragment — which is often 
highly negatively charged — have a tendency to 
extend far out of the fragment volume, in region 
of space where the other electrons of the crystal 
should be located. It is therefore necessary to for- 
bid these domains of 1Z 3 by a repulsive potential 
mimicking the electronic volume of the atoms of 
the rest of the crystal. This aim can be reached by 
the usage of total ions pseudo-potentials 21 at the 
location of the first shells neighboring atoms of the 
fragment. 

Since all these embedding effects are mono- 
electronic they can be easily checked. Indeed, they 
are included in a simple Hartree-Fock (HF) calcu- 
lation thus a good embedding should reproduce the 
HF density of states of the entire crystal, projected 
onto the fragment atomic orbitals. 
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4-A. Example : the a'NaV^Os compound 

As an example of the importance of the embed- 
ding on the orbitals definitions and orbitals en- 
ergies 22 we will discuss the case of the a'Nal^O^ 
compound. a'NaX^O^ is a ionic layered compound 
built from V0 5 square-pyramids arranged in (a, b) 
planes (see figure EJ) and doped by sodium ions 
located between the layers. 






FIG. 9: Schematic structure of a AfaViOs, along the a-b and 
the a-c planes. The oxygen atoms are denoted by open circles, 
the vanadium atoms by filled circles and the sodium atoms by 
dots. 

The compound is strongly ionic and the 
formal charge transfer yields a formula of 
Na + (V 2 ) 9+ (0 2 ~) 5 . The vanadium atoms support 
therefore one unpaired electron for two magnetic 
sites. It has been shown that the different layers 
are very weakly coupled and that the electronic 
structure can be considered as two dimensional 2 ^. 
In fact most authors see the compound as coupled, 
quarter-filled (one active electron for two active or- 
bitals), magnetic ladders (see figure ITUj) . 

Let us concentrate, for instance, on the interac- 
tions of the magnetic electrons on a rung of the 
ladders (for an ab initio evaluation of the differ- 
ent interactions in the high temperature phase see 
ref.—, in the low temperature phase see ref.— and 
for an ab-initio evaluation of the charge ordering 
ref.—). They are of two types, hopping and super- 
exchange. They can be evaluated by the spec- 



FIG. 10: Schematic structure of the coupled magnetic ladders 
in the (a,b) plane. 



troscopy of an embedded fragment composed of 
two VO5 pyramids sharing a corner oxygen (see 
figures HU and H2J) . 




FIG. 11: Fragment of the crystal used for the cluster calculations 
of the interactions along a rung seen in the (a, b) plane. 




FIG. 12: Fragment of the crystal used for the cluster calculations 
of the interactions along a rung seen in the (a, c) plane. 



We can therefore check the importance of the 
embedding on the orbitals definitions and orbitals 
energies by the comparison of the density of states 
(DOS) of the entire crystal (as computed from pe- 
riodic HF techniques 27 ), the density of states of 
the isolate fragment and the density of states of 
the fragment embedded in different baths. For 
sake of clarity we will project the DOS on the 
most important orbitals for the effective hopping 
and exchange interactions, namely the magnetic 
orbitals of the vanadium atoms (3d xy if x, y, z re- 
fer respectively to the crystallographic orthogonal 
a, b, c directions) and the 2p y orbital of the bridging 
oxygen that mediate the interactions. Figure [T3] 
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and figure ITU report respectively the different pro- 
jected DOS on the magnetic and bridging orbitals. 
One can see immediately that the isolated frag- 
ment DOS is very different from the infinite crystal 
one. In fact, while the highest occupied crystalline 
bands are built from the d xy atomic orbitals of the 
vanadium atoms, in accordance with the magnetic 
character of this orbital, in the isolated fragment 
the molecular orbital based on the vanadium 3d xy 
is far above the Fermi level, the unpaired Highest 
Occupied Molecular orbital being a delocalized or- 
bital on the 3p orbitals of the vanadium and the 2p 
orbitals of different oxygen atoms. The embedding 
built from atomic formal charges (V 4,5+ , 2 ~ and 
Na + ) is properly reproducing the DOS projected 
on the oxygen 2p y bridging orbital, however it is 
unable to reproduce the crystal DOS projected on 
the magnetic orbitals. An orbital analysis of the 
different embedded fragments and crystal calcula- 
tions shows that the apical oxygen is not linked to 
the vanadium by an ionic bond but that a strong, 
dative, triple bond is taking place between these 
two atoms. The consequence is a strong devia- 
tion of the Madelung potential from the one de- 
fined by formal charges. A redefinition of the bath 
using approximate self-consistent Miilliken charges 
on the vanadyle, that is V 3+ and O 5 ~ finally yield 
a proper reproduction of the DOS projected on all 
the orbitals of the fragment, insuring a reliable tak- 
ing into account of the non-local crystal effects. 



Crystal 



Ai. 



Mulliken charges 



Formal charges 



-0.6 -0.2 0.2 

FIG. 13: Density of states 
projected on the magnetic 
d X y orbitals of the rung 
vanadium atoma— . 



Crystal 



Mulliken charges 



Formal charges 




-0.6 -0.2 0.2 



FIG. 14: Density of states 
projected on the 2p y orbital 
of the rung bridging oxygen 

, 22 

atorr— 



Let us now take a look at the importance of 
these embedding problems on the effective inte- 
grals amplitudes, namely the hopping t± integral 
of a magnetic electron between the two vanadium 



atoms of the rung, and the super-exchange J± in- 
tegral of two magnetic electrons on the rung. As 
seen previously, these integrals can be defined from 
the fragment spectroscopy as half the symmetric 
versus antisymmetric doublet excitation energy for 
the hopping (t ± = 1/2 \E(D g ) - E(D U )]), and as 
the singlet versus triplet excitation energy for the 
exchange (J± = E{Sg) —E(Tp)). We shall not de- 
scribe here the computational technique used for 
these spectroscopic calculations since the next sec- 
tion will be devoted to this problem. Table 1 inte- 
grals both for the formal charges and the Madelung 
charges embeddings. One can see immediately 
that a small change in the dipole moment of the 
vanadyle fragment (V 4 ' 5+ 2 ~ in the formal charges 
bath, V 3+ 0°' 5 ~ in the so-called Mulliken charges 
bath) can considerably modify the effective inte- 
grals amplitudes, since J± is multiplied by a 1.73 
factor between the two calculations. 



Integral 


Mulliken charges 


Formal charges 


J± 


-293.5 meV 


-509.7 meV 


t± 


-538.2 meV 


-581.2 meV 



TABLE 1. Hopping and super-exchange effective integrals 
of magnetic electrons between the two vanadium sites of 
the rung according to the nature of the embedding^. 

5. Fragment spectroscopy 

The fragment embedding being properly defined, 
the question is now how to computed accurate ex- 
citation energies and related wave-functions. The 
question is far from trivial since the total energies 
of the computed states are usually of the order of 
10 3 eV A while the seeked transition energies are of 
the order of a few tenth up to a few tens of meV, 
that is from 5 to 7 orders of magnitude smaller. 
Even the correlation part of the energies are of the 
order of a few tens of eV which is still between 
3 to 5 orders of magnitude larger than the seeked 
excitations energies. Reliable results therefore ne- 
cessitate procedures that calculate identically the 
major part of the total energies. This can only be 
achieved if a common set of unbiased orbitals is 
used for the different states involved in the exci- 
tation, whether they do belong to the same irre- 
ducible symmetry representation (space and spin) 
or they do not. 

A zeroth-order description of these seeked tar- 
get states can be obtained using a CASCI 38 cal- 
culation. The main question is the choice of the 
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set of orbitals supporting the CAS space. While 
it is obvious that the active orbitals (as defined 
in sections 3 and 1) should belong to this set, 
one may wonder whether it should be restricted 
to them. The question is specially crucial in the 
case of bridged systems where the interactions are 
mediated by a ligand and the question of includ- 
ing in the active space the mediating orbitals of the 
bridge (denoted (IS X ) should be asked. The answer 
depends on the weight — in the target states - 
of the configurations involving a partial occupa- 
tion of x ; that is different from 2 if x is a core 
orbital and different from if virtual one. If 

large, the mediating orbitals x should be included 
into the set supporting the CAS. In such cases, one 
may wonder whether the definition of the active or- 
bitals supporting the model shouldn't be enlarged 
accordingly. Independently of the definition of the 
CAS for the ab-initio calculations, the answer to 
this question and accordingly the validity of the 
model and related active orbital space definitions, 
can be checked using a simple criterium : the norm 
of the projection onto the active space supporting 
the model Hamiltonian, of the computed target 
states wave functions. A large value (> 0.7) of 
this projection insures a correct definition of the 
model while a small value is the sign that an im- 
portant physical effect have been omitted. 

While the CAS is a good reference description, 
it is well known that a CASCI calculation is far 
from yielding accurate excitation energies. Indeed 
it is missing not only the dynamical correlation ef- 
fects, but also the space and spin repolarisation 
effects specific to the different VB configurations 
(dynamical polarisation), which are known to be 
crucial both for a proper treatment of charge delo- 
calisation 28 and for correct relative weights of the 
different VB configurations (in particular the ionic 
ones) in the active spaced- . This last point is of 
great interest since it discards all computational 
methods which do not allow a revision of the rel- 
ative weights of the different VB configurations 
within the CAS, in perspective of the dynamical 
polarisation effects. This is the case for instance 
of the perturbative methods after a CAS zeroth 
order. They necessitate large CAS including the 
major part of the repolarisation effects in the ze- 
roth order wave function in order to give accurate 
results^. 

As a matter of illustration table 2 reports the 
evaluation of the effective exchange and hopping 
integrals on the rung of the a'Nal^O^ compound 
computed both at the CASCI (the CAS being 



supported by the magnetic orbitals and the me- 
diating orbital of the bridging oxygen) and the 
CASCI + all single excitations on all CAS deter- 
minants (zeroth-order + dynamical polarisation) 
levels. One notices immediately that the effective 
integrals are dramatically modified by the inclu- 
sion of the dynamical polarisation effects. Indeed 
their inclusion increases the super-exchange inte- 
gral by a factor ~ 5.3 and the effective hopping 
by a factor of ~ 1.3. One may however wonder 
why this repolarisation effect is so much stronger 
on the exchange integral than on the hopping. At 
the simplest level of description, the hopping in- 
volves a charge fluctuation between the two mag- 
netic sites, thus strong changes in the electrostatic 
potential. An accordingly strong dynamical po- 
larisation response can be expected. On the con- 
trary the exchange integrals involves only spin field 
modifications which is usually expected to yields 
much smaller repolarisation effects. This apparent 
contradiction can be levelled if one remembers the 
through-bridge super-exchange mechanism shown 
in figure 121 While the through-bridge hopping 
mechanism involves only paths going through con- 
figurations with a le occupation change (relative 
to the references) of the bridging orbital x, the ex- 
change mechanism involves not only similar paths, 
but also paths where the bridging orbital occupa- 
tion is modified by 2e. This is these paths that 
induce extremely strong specific repolarisation ef- 
fects. Indeed, in the present case this last path 
carries half of the effective exchange amplitude ex- 
plaining the observed much stronger repolarisation 
effects than on the hopping integral. 



Integrals 


CASCI 


CASCI 


CASCI 






+ singles 


+ singles 








+ diff. doubles 


J± 


-60.6 


-321.1 


-293.5 


t± 


-420.6 


-542.7 


-538.2 



TABLE 2. a'NaV^Os. Effective exchange and hopping 
integrals between the magnetic d xy orbitals of the two vana- 
dium atoms of the rungs24. The calculations have been 
performed using the IDDCI„ approach. 

Comparatively the dynamical correlation effects 
due to the double-excitations on the CAS determi- 
nants have much smaller contributions (in general 
of the order of 10% to 20%2i and much smaller 
in the present example). I would like to point 
out that while the large majority of the double- 
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excitations on the CAS involve excitations from 
two core orbitals towards two virtual orbitals, 
these excitations do not significantly contribute 
to the excitation energies between two CAS-based 
states. In fact they exactly cancel out in a sec- 
ond order perturbative approach. It can therefore 
be a wise solution to consider only the differen- 
tial double-excitations, that is only the double- 
excitations that contribute to the energy differ- 
ences at the second order of perturbation. This 
is what is done in the Difference Dedicated Con- 
figuration Interactions method of Malrieu et a@&. 

One should however notice that in the presented 
results the bridging orbital have been included in 
the CAS and that the crucial bridge polarisation 
and correlation effects are included in the single 
and double-excitations to the CAS. Whether the 
CAS should be restricted to the only magnetic or- 
bitals, than triple-excitations to the CAS should 
be considered in order to treat these effects. 

6. From the ab initio wave 
function to the model 
Hamiltonian 

In the simplest cases — which are the most nu- 
merous — the knowledge of the excitation energies 
is necessary and sufficient to determine the inte- 
grals amplitudes of the model hamiltonian. This 
is for instance the case for the effective exchange 
integral between 1/2 effective spins (as was largely 
discussed in section 2) or for the effective hopping 
between sites, if they are not equivalent by the 
space symmetry operations. However things can 
be sometimes more complicated. Typically one 
may encounter two types of difficulties. 

• The problem may be over-determined. 
This is for instance the case for an effective 
spin-1 Heisenberg model. Such models are 
particularly adequate for the description of 
magnetic systems where the magnetic sites 
support two unpaired electrons in (quasi-) de- 
generate, strongly correlated, orbitals. The 
Hund's rule imposes a triplet coupling of these 
two electrons and, under certain conditions^, 
they can be modeled by an effective spin- 
1. The interaction between two such spin- 
1 yields a set of singlet, triplet and quintet 
states, which wave functions are totally de- 
fined by the spin symmetry, and energy dif- 
ferences by the effective exchange parameter 



J. One therefore has two computed target 
states excitation energies (singlet to triplet, 
AE S g^t P , and singlet to quintet, AE sg _, qt ) 
that have no reasons to strictly follow the 
spin-1 Heisenberg hamiltonian eigenvalues re- 
lationships (AE sg ^ tp = J, AE sg _ qt = 3 J), 
and only one parameter (J). 

• The problem may be under-determined. 
This is for instance the case for a hopping 
problem between non-equivalent sites : the 
excitation energy between the two doublet 
states is not sufficient to determine both the 
inter-site hopping and the energy difference 
between the two active orbitals. In such cases, 
it is necessary to use the wave functions in- 
formation. 

• More generally in systems for which the states 
are not totally determined by the space or 
spin symmetry, the energy differences are not 
sufficient for the complete and accurate deter- 
mination of the effective model. Since the ex- 
citation energies cannot provide the necessary 
wave function information, it is impossible to 
meet the requirement of a good reproduction 
of the wave-functions physical content. 

According to the types of difficulties encountered, 
several solutions can be proposed that will be de- 
tailed in the next sections. 

6-A. Least square fit 

In the case of an over-determined problem, one 
may choose to abandon the reproduction of some 
target states and to reduce the objectives of the 
model to the accurate description of a smaller set 
of selected states. Part of the model hamiltonian 
eigenstates will thus be present in the model only 
as a purpose of completion or in order to allow the 
sufficient flexibility to reach the seeked accuracy 
on the selected states. A detailed analysis of this 
option will be given later on. 

An alternative choice can be to abandon the ex- 
act reproduction of the ab initio data and rather 
to model them at the best using a least square fit 
technique. This method has the advantages to de- 
scribe equally well (or bad) the whole set of tar- 
get states. The aim is to fit, simultaneously and 
on an equal footing, the computed energies and 
target states projections. One will thus apply the 
least square fit on the eigen-equations of the model 
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hamiltonian (Hmd). In a first approximation the 
Lagrangian can be denned as 

i 

or equivalently 

£ = EE(( $ «l^ p l*i>/(*«l p l**>- e i) 2 

i a 

where |^) and £{ are respectively the computed 
target states wave functions and energies, P = 
J2 a \ ®a){$a\ is the active space projector and 
are the VB configurations acting as basis set of the 
active space. The minimisation of £ runs on the 
parameters of Hmai- The strict usage of the above 
formulation presents however some drawbacks. 

• The total energies have no real signification 
and only the excitation energies have a phys- 
ical meaning. One should therefore use only 
eigen-equations differences in the Lagrangian. 

• The weight devoted in the Lagrangian to the 
very small components of the wave function 
is much too large compared to the weight 
devoted to the large components. Specially 
since the relative error made on these small 
coefficient is much larger and their physical 
importance much less. One should there- 
fore ponderate their contribution in L so that 
small variations of their amplitude (but that 
may be large in relative variations) do not sig- 
nificantly affect the result. 

Finally, as in all least square fit method, the accu- 
racy of the procedure is measured by amplitude of 
\/C which should be much smaller than the typical 
scale of the fitted effective interactions. H-md 

6-B. Using the wave function information 

The usage of the computed target states wave 
functions raises the question of the reduction of the 
information contained in a wave function spanned 
over several hundred thousands to millions of de- 
terminants onto a model state described by a few 
VB configurations. 

For this purpose it is possible to use the guide- 
lines of the effective hamiltonian theory^, that is 
to define the model hamiltonian by its eigenvalues 
and eigenvectors. 

• The model hamiltonian eigenvalues (e$) 
should be set to the target states excitation 
energies computed from the ab-initio spec- 
troscopy of the embedded fragments (eA 



• The model hamiltonian eigenvectors 
should be set to the projection onto the ac- 
tive space of the computed target states wave 
functions 

These conditions — that suppose the equality be- 
tween the active space and the target space dimen- 
sions (m a = m t ) — totally define the model hamil- 
tonian in a unique manner. However the procedure 
can be strictly used only in few special cases since 
it is often incompatible with some the of formal 
requirements demanded from the model. 

• Typically when the target states projections 
onto the active space (P|\&j)) are not fully de- 
fined by the symmetry of the system, there is 
no reason for them to form an orthogonal set. 
The direct corollary of this non-orthogonality 
is the lost of the hermiticity property of the 
model hamiltonian. Going back to the exam- 
ple of the effective hopping between two non 
equivalent sites a and b, the two computed 
target doublets (|^i) of energy e\ and of 
energy £2) can be written as 



P|*!> 



cos a I a) + sin a \b) 
— sin (3 I a) + cos (3 \b) 



which yield a model hamiltonian 



/ £1 cos a cos (3+E2 sin a sin ft 
" cos (a— ft) 

(ei —£2) sin a cos ft 
\ cos (a— ft) 



\b) 

(ei — £2) cos a sin/3 
cos (a— ft) 

£1 sin a sin ft+62 cos a cos ft 
cos (a— ft) 



One sees immediately that when a 7^ (3 the 
hermiticity is not verified. 

• In addition, the above procedure yields a pri- 
ori independent and non- zero matrix elements 
between all VB configurations in the same ir- 
reducible representation. In consequence, the 
number of independent parameters of the ef- 
fective model is equal to the number of matrix 
elements not directly imposed by the symme- 
try. The goal of the modelisation, in terms of 
the reduction of the physics to a small number 
of physically readable, dominant interactions, 
is therefore totally lost. 

The price to retrieve the hermiticity and the sim- 
plicity is either (i) the abandonment of the exact 
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reproduction of the computed excitation energies 
and wave functions to the profit of an approximate 
representation through a least square fit method 
(see previous section), or (ii) the abandonment of 
the reproduction of part of the target states in or- 
der to focus on the accurate description of a set of 
selected states. In the above example of asymmet- 
ric hopping the choices would translate into either 
the approximate representation of the two doublets 
by average states as resulting from a Lowdin or- 
thogonalisation procedure 



P\*2) 

(WO 



cos — - — \a) + sin — - — \o) 



— sin 



2 

a + /3 



\ _l a + f 3 \h\ 
a) + cos \b) 



or the accurate reproduction of one of the doublets 
(for instance |^i)), the second one being accord- 
ingly defined for this purpose (Smidt orthogonali- 
sation) 



(*i\P\*i) 

(*l\P\*l) 



cos ol \a) + sin a \b) 
— sin a | a) + cos a \ b) 



More generally this second choice is related to 
the guidelines of the intermediate hamiltonian the- 
ory^ : reproduction by a subset of the active 
space of a selected number of target states, both in 
terms of energies and wave functions. One should 
therefore partition the active space underlying the 
model hamiltonian into two subspaces. 

• A main active space of dimension (m ma < 
m a ). This subspace will be supported by the 
eigenstates of the model hamiltonian aimed 
at reproducing the target states. m ma should 
therefore be equal to the number of target 
states to reproduce {m ma = m t ). 

• An intermediate active space of dimension 
m ia = m a —m ma , supported by the remaining 
eigenstates of the model hamiltonian. These 
states are devoted to provide to the model the 
necessary flexibility to both verify the (for- 
mal) requirements and to best reproduce the 
target states. They have a formal purpose but 
no physical purpose even though they may be 
related — however not a trustworthy repre- 
sentation — to meaningful physical states. 



Such a procedure, based on a selective choice, send 
us back to the aim, in fine, of the model. Indeed 
the partition of the active space must be deter- 
mined by the physics targeted at, that is the states 
of the entire lattice one would like to describe. 
From the embedded fragment point of view, it 
means an accurate reproduction of the local phys- 
ical content of these infinite system states, that 
is of the fragment states which are highly popu- 
lated (as can be computed by the density matrices 
projected on the fragment) in the infinite system 
states. 

Let us illustrate our purpose by a simple exam- 
ple : the low energy physics of a half-filled mag- 
netic system. The infinite system is therefore a 
lattice which low energy physics is based on equiv- 
alent magnetic centers supporting each one active 
orbital and one active electron. The states respon- 
sible for the magnetic properties are therefore neu- 
tral in a Valence-Bond sense, that is configurations 
where each magnetic center support a unique ac- 
tive electron are dominant in the low energy states 
of the entire system. Considering now an embed- 
ded fragment containing two magnetic centers (say 
A and B), the active orbitals (a and b) and elec- 
trons define a 2e, 2 orbitals CAS underlying a Hub- 
bard model hamiltonian. 
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of energy U 



In a Bloch effective hamiltonian theory the above 
problem is over-determined since one has 3 energy 
differences and 2 singlet eigenstates which wave 
functions are not totally determined by the sym- 
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metry. The result is 5 conditions for only 2 un- 
known parameters. The infinite system physics 
being supported by neutral VB configurations it 
is crucial that the local neutral states are best re- 
produced. One should therefore prefer an inter- 
mediate hamiltonian guideline and focus on the 
reproduction of the \Sg) and \Tp) energy differ- 
ence and wave functions, yielding a well defined 
problem with two conditions (reproduction of the 
singlet-triplet energy difference and ionic config- 
urations coefficient in the singlet state) and two 
unknowns. The main consequence of this choice 
is that the on-site effective bi-electronic repulsion 
integral U is determined so that to best reproduce 
the correct component in the singlet ground state 
of the ionic VB configuration \aa) + \bb). As a 
corollary the energies of the essentially ionic states 
\Su) and \Sg*) are overestimated. Indeed, it can 
be seen from a simple perturbative expansion that 
the dynamical correlation effects on the ionic con- 
figurations strongly depend on whether they be- 
long to an essentially neutral or essentially ionic 
state. In case of a neutral state the rescaling of 
the ionic configurations energy by the dynamical 
correlation processes is (at the second order of per- 
turbation) 



U 



\(aa\l/r 12 \^)\' 
Ei 



while in case of an ionic state it comes 



U 



\{w\l /r 12 \$j 
Ei-U 



One should thus be conscious that while the model 
determined under the previous criteria will be ap- 
propriate for the determination of the low energy 
physics it is unsuitable to study the transitions to- 
ward essentially ionic VB states. 

In this example the reduction of the number 
of target states to reproduce provides a solution 
to two different problems. The first one is get- 
ting a mathematically well conditioned inversion 
problem. The second one is the resolution of the 
dilemma between an accurate reproduction of the 
ground state wave function and of the transition 
energy towards ionic excited states. While the sec- 
ond type of problems can be solved by the usage 
of intermediate hamiltonian guidelines, the solu- 
tion to first type of problems may require such 
drastic reduction of the target space (one state 
per irreducible representation) that the model may 
become meaningless. In such thoughtful 



combination with square fit techniques will provide 
the exit. 

7. Conclusion 

The fascinating properties observed in the new 
materials synthesized in the last couple of decades 
are commanded by their electronic structure and 
more specifically by strong electron-electron cor- 
relations. The understanding of the observed col- 
lective effects necessitates a clear picture of the 
interplay between the dominant microscopic local 
interactions. Quantum chemists with their exten- 
sive culture of correlation effects in finite systems, 
of valence correlation but also dynamical correla- 
tion and its retroactive effects on the valence shell, 
in the ground state but also in excited, multi- 
reference states, have the intellectual tools to an- 
swer this question. We have seen in this paper how 
up-to-date quantum chemical spectroscopic meth- 
ods associated to a controlled description of the 
crystalline embedding, a clear understanding of 
the model hamiltonians aims, an appropriate ex- 
tension and usage of effective hamiltonian methods 
that allow a controlled reduction of the ab initio 
information into a few dominant interactions, are 
the clues to a controlled and accurate modelisation 
of these materials. 

It is a fascinating challenge for the quantum 
chemist used to small finite systems, to enter this 
condensed matter field where the answers can only 
be find by the integration over all physical scales 
from first principle microscopic local electronic in- 
teractions, to the nanoscopic scale (fragments) up 
to the macroscopic properties. 
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The extended Hubbard model 



The Heisenberg model 



The Heisenberg model was first devised by 
Heisenbergii and formalized by Van Vleck and 
Dirae— in order to rationalize isotropic ferro- 
magnetic interactions between localized electrons 
or spin 1/2 particles supported by orbitals lo- 
calized on well separated centers. Its validity 
has later been extended by Anderson^ to anti- 
ferromagnetic interactions by the famous super- 
exchange mechanism : stabilisation of the singlet 
state compared to the triplet state by the effec- 
tive effect of the ionic configurations. The general 
expression of the Heisenberg Hamiltonian is for a 
system with one spin, one orbital per center (half- 
filling) 

H He is = Jij (Si-Sj - jniUj] (1) 
<ij> V 4 / 

When the spins represent electrons (spin 1/2) the 
Heisenberg Hamiltonian can be rewritten as 

= j ( a \Al- a )Al \ ( gjt%l ~ gjTM oh 
eiS <tT> 13 \ v/2 ) \ V2 7 

where a\ a (resp. a ia ) is the creation (resp. annihi- 
lation) operator of an electron of spin a on site i, 
rii is the number operator on site % and S, the spin 
operator on site i. The sum over < i, j > runs on 
nearest neighbors sites i and j. 



The extended Hubbard model is basically a sim- 
plified version of the Pariser-Parr-Pople hamilto- 
nian. It contains essentially four types of inter- 
action parameters, two mono-electronic ones cor- 
responding to the Hiickel approximation and two 
bielectronic ones 

• the orbitals energies, 

• the hopping integrals (usually limited to ad- 
jacent centers), 

• the bielectronic coulomb repulsion between 
two electrons located on the same orbital, 

• and the bielectronic coulomb repulsion be- 
tween two electrons located on different or- 
bitals (most of the time limited to nearest 
neighbors centers). 

The three first terms define the Hubbard model 
and the addition of the inter-site repulsions leads 
to the extended Hubbard model. It comes the fol- 
lowing generic formulation 

H Hub — X! 5Z ( a !o- a i°- + a ]a a icr) + 

<ij> & 

Y UiH^nn + (4) 

i 

Y Vij (m - Zi) (nj - Zj) 

where Zj is the effective nucleus charge of the cen- 
ter i. 



The t — J model 

The t — J model was first devised by Zhang and 
Ricei^ in order to described the movement of holes 
in an anti-ferromagnetic background such as found 
in the copper oxide planes of the cuprate high T c 
super-conductors. It applies to magnetic systems 
away from half-filling (number of magnetic orbitals 
different from the number of electrons or spins). 
According to the sign of the doping (extra elec- 
trons or extra holes compared to the half-filled 
system) the hamiltonian applies to the electrons 
or holes. Its general form is 

Ht-j = E J H ( s *- s i - J n i n j) + 

<i,j> V ' 

J2 hi J2 ( a t%<x + a) a a ia j (3) 

<i,j> ° 
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